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ABSTRACT 

The new photometric space-borne survey missions CoRoT and Kepler will be 
able to detect minute flux variations in binary stars due to relativistic beaming 
caused by the line-of-sight motion of their components. In all but very short 
period binaries (P > 10 d), these variations will dominate over the ellipsoidal 
and reflection periodic variability. Thus, CoRoT and Kepler will discover a new 
observational class: photometric beaming binary stars. We examine this new 
category and the information that the photometric variations can provide. The 
variations that result from the observatory heliocentric velocity can be used to 
extract some spectral information even for single stars. 

Subject headings: techniques: photometric — techniques: radial velocities - 
surveys — binaries: close — binaries: general — binaries: spectroscopic 

1. Introduction 



In 2003 Loeb fc Gaudil suggested a new photometric method to detect extrasolar planets, 



based on the minute variab ility of the stellar flux due t o relativistic beaming induced by the 



star's reflex radial velocity. iRybicki fc Lightmanl (119791 ) show that several factors contribute 



to the beaming effect. The bolometric factors are the Lorentz transformations of the radiated 
energy and the time intervals, and the modification of the angular distribution of the radiated 
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energy (stellar aberration). Thus, in the limit where the star's radial (line-of-sight) velocity 
v K is much smaller than the speed of light c, the observed bolometric flux F is modified 
relative to the emitted bolometric flux F as 



1+4— 

c 



X 



For band-pass photometry, the Doppler shift of the emitted frequency and the spectral index 
of the source spectrum also have to be taken into account, which finally yields: 



F v 



(3 — a) — 

c 



(2) 



where F v and F u0 are the observed and emitted flux density at frequency u, respectively, and 
a = d log F/d log v is the average spectral index around the observed frequency. Note that 
Equations [1] and [2] can also be obtained in a semi-classical context. In any case, since the 
effect depends on the first order of v R /c the velocities do not need to be highly relativistic 
to detect the effect. Detecting beaming variability related to orbital motion requires the 
very high precision provided by photometric satellites. Even in the extremely high orbital 
velocities of the stars orbiting the massive black hole in the Galactic Center, the beaming 
effects cannot be m easured with the low photometric precision available for those stars 
rtZucker et aDl2006h . 



Loeb fc Gaudil (120031 ) suggested to detect stellar periodic motion induced by an unseen 



planet through the beaming effect. They showed that the amplitude of the beaming ef- 
fect produced by an extrasolar planet is of the order of micro- mag nitude . Such an effect 



is ba rely detectable by the new ph otometric space missions CoRoT (IBagliru I2003T ) and Ke- 
pler (IBasri. Borucki fc Koch! 120051 ). which are aiming to find extrasolar tra nsiting planets, 
with a typical variability amplitude of 100/zmag. The obvious advantage of lLoeb fc Gaudi 
( 120031 ) approach is its applicability to planets with almost any inclination, whereas detecting 
transiting planets is limited only to planets with orbital inclinations close to 90°. 

Here we examine the beaming effect in binary stars, and suggest a new class of binaries 
- beaming binaries. These are a hybrid between spectroscopic and ellipsoidal binaries, since 
the beaming binaries will be detected by periodic photometric variations due to their orbital 
radial velocity. We show that a binary with an orbital period of P = 100 d has a beaming 
variability of at least (9(100 fimag), which is easily detectable by CoRoT and even more so 
by Kepler. We further show that for binaries with periods longer than 10 days, beaming 
variability dominates over the competing effects of ellipsoidal and reflection variability. We 
therefore expect the new satellites to harvest hundreds of previously unknown binaries of 
this new class. 



Beaming variability of order 100 //mag will also be induced by the Earth's heliocentric 
motion relative to any observed source, single star or binary. In spectroscopic observations, 
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the effect of the Earth motion is corrected in order to obtain heliocentric radial velocities. 
Here we show how the dependence of this small variability on the spectral slope at the 
observed bandpass can be used to probe the spectral characteristics of single and binary 
stars. 



2. The amplitude of the beaming effect as compared with the ellipsoidal and 

reflection effects 

The three kinds of periodic flux variations we expect to detect in binary stars are those 
due to the ellipsoidal tidal deformations of the stars, the reflection of the lig ht of each star 



by its companion, and relativistic beaming, which is the subject of this work. lLoeb fc Gaudi 



( 120031 ) have presented a rough comparison of those three effects in the context of planet- 
hosting stars. We now compare the three effects for binary stars, assuming for simplicity 
a circular (e = 0) edge-on (sin 2 = 90°) orbit. Note that for systems with two stellar 
components, the observed ellipsoidal effect is the weighted average of the ellipsoidal effects 
of the two stars, whereas the observed beaming effect is the weighted difference between the 
beaming variabilities of the two stars, as their effects are exactly in opposite phase. In this 
respect beaming is similar to the reflection effect. 

To a good approximation (see below), the magnitude of the beaming effect can be 
calculated under the assumption that the two stars radiate as blackbodies. For a blackbody 
source of temperature T e g the spectral index a is: 

e x — 1 

where x = hv/kT e g. 

The binary's orbital separation is given by Kepler's third law, 



Mi + M 2 \ 1/3 / P \ 2/3 



M J VlyearJ 



AU , (4) 



where M, are the two masses, with the subscript 1 referring to the primary and 2 to the 
secondary. The amplitude of the primary's radial-velocity variation is then 

Kl= ( T7 Ml Tr )( ' r 1/2 29.8hns-' . (5) 



M 1 + M 2 J VI AU 



A corresponding expression is obtained for the secondary by interchanging the subscripts 1 
and 2. The peak-to-peak amplitude of the total expected relative flux variation from the 
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binary due to beaming is then 

fAF u \ _ 1 ^[3 - ai^)]^ - K 2 [3 - a 2 (^)]F z 



V F J c F +F (6) 

\ 1 v / beaming L " ^ f,l T -£ „ )2 

Note that because the observed effect is the difference between the effects of the two stars, 
the beaming effect vanishes for an equal-mass binary, where the spectral characteristics are 
also identical for the two components. 

In order to e s timat e the ellipsoidal variability, we use the expression presented by 



Morris fc Naftilanl (119931 ) for the peak-to-peak ellipsoidal variability of the primary: 



AF„A (15 + m)(l + T0 M a (RiV 

^Jellips ' 3-Ul MiU/ 

Here T\ is the gravity darkening coefficient of the primary and u\ is its limb-darkening 
coefficient. We calculate a similar expression for the secondary and then weight them by the 
expected blackbody fluxes in order to obtain the total relative variation of the binary. 



Morris fc Naftilanl (119931 ) also provide a prescription for calculating the amplitude of 
the reflection effect. They assume that each star absorbs some of the bolometric flux of 
its companion, which heats the stellar hemisphere facing the co mpanion, inducing a n asym- 



metric increased emission. Assuming blackbody radiation law, iMorris &: Naftilanl define a 
'luminous-efficiency' factor by: 

, ,'T 2 \V 2 -1 



T-i 



XI _ 1 ' 



where T\ and T 2 are the temperatures of the two components. Like the beaming effect, the 
contributions of the reflection effects of the two stars are in opposite phase, and the total 
magnitude of the effect is the weighted difference of the two. Keeping only the leading order 
terms in the radii (expressed in terms of the orbital separation), we obtain: 

AF U \ _2(R 2 /aff- 1 F u , 1 - (Rja) 2 f x F u>2 



J reflect ^ ^V,l + -^,2 



Note that the ellipsoidal and reflection variabilities were calculated assuming tidal lock- 
ing and a circular orbit. Furthermore, we use only the leading order terms in the fractional 
radii. Thus, we might be overestimating the amplitude of those effects. However, these 
expressions suffice as conservative estimates for comparing with the beaming effects. 

We use Equations [61 dand [9] to compare the three effects for three typical binaries and 
for a range of periods. Table [1] presents the parameters assumed for the stellar components. 
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For the purpose of this simple comparison, we assume U\, 112 = 0.6 for ou r hypothetical stars 



We calculate the gravity darkening coefficients using the prescription in iMorrisI (119851 ) 



Table [2] compares the three effects for the three binaries, observed in the V band, for 
periods of 10 and 100 days. In all cases the beaming variability dominates over the other 
two effects. Figure [1] shows the three effects for an F0-K0 binary for a range of periods. The 
corresponding plots for the other cases were very similar. The dependence of the effects on 
the orbital separation is explicit in the expressions above, and we can use it to understand 
the dependence on the orbital period. While ellipsoidal variability decreases with period as 
P~ 2 , and the reflection variability as P~ 4//3 , the beaming variability only decreases as P -1 / 3 , 
and we expect it to become dominant for long enough periods. 

In the three cases we examined the ellipsoidal variability dominates for periods shorter 
than 8 days, while for periods longer than 10 days the beaming variability becomes dominant. 
Remarkably, the three lines intersect at about the same period, and the reflection effect is 
almost never dominant. Furthermore, the amplitude of the beaming variability stays at the 
detectable levels for C0R0T and Kepler, of 0.1 — lmmag for periods of a hundred days and 
more. 

In Figure [2] we compare the three effects for a range of mass ratios. We assumed a GO 
primary, and used power laws for the dependence of the secondary radius and temperature on 
its mass: R oc M a8 and T e g cx M 55 . While the ellipsoidal variability is mostly sensitive to 
the highest mass ratios, the beaming and reflection effects are more sensitive to intermediary 
mass ratios. This is mainly because at the highest mass ratios the effects from both binary 
components are canceled out. 



3. Discussion 

The radial-velocity beaming lightcurve can yield directly most of the spectroscopic or- 
bital elements, including the period, eccentricity and time of periastron passage. Since these 
will be obtained as the result of a well defined magnitude-limited photometric survey, they 
will provide large amounts of new data to statistical studies of spectroscopic binaries, includ- 



ing, e .g., the distribution of orbital period (jDuqennov fc Mavorlll99ll;lMazeh. Tamuz fc North 



20061 ). and the relation between orbital period and eccentricity (jHalbwachs et al.l 120031 ). 



The only spectroscopic element that cannot be obtained directly from the lightcurve 
is the radial- velocity amplitude K\. However, as outlined in Equation [6], the K\ value can 
be derived from the amplitude of the beaming effect through the spectral index a of the 
primary and the relative amplitudes of the beaming effect of the two components of the 
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binary. In most binaries, the secondary is faint enough that we will be able to ascribe the 
observed beaming variability solely to the primary component. If the primary spectral type 
is known, we can derive its spectral index (a) and obtain K lt thus deriving the full set of 
orbital elements of a single-lined spectroscopic binary. 

In order to estimate a and calibrate the relation between the beaming amplitude and 
K\, some spectroscopic follow-up observations of the detected binaries should be performed. 
Since most of the radial-velocity elements will already be known from photometry, only 
a small number of observations is ne eded per star. Mult i-object spectrographs, such as 



FLAMES on the Very Large T elescope (jPasquini et al.ll2002l ) or Hydra on the Wiyn telescope 



( IBarden fc Armandrofflll995l ). seem like an efficient means to obtain these observations for 



the detected beaming binaries in the field. 

In the few cases where the two components might have very similar magnitudes and 
masses, the two contributions to the beaming variability may cancel out, because of their 
opposite phases. In cases where the secondary light will be significant but will not cancel 
the primary light completely, we will need a photometric analogue of spectroscopic disentan- 



gling procedures like TODCOR (IZucker fc Mazehl Il994l ). Measurements in more than one 



photometric band may add the constraints needed to solve for K\ and K-i- 

We note in passing, that the derivation of K%, and when possible also K2, is sensitive to 
any blending of the binary image with other stars, as they depend on the relative amplitude 
of the beaming effect. Therefore, it would be needed to obtain high resolution image of 
the observed field, in order to spot any other possible contributions to the binary light that 
might dilute the beaming effect. In fact, measurements in different bands may also serve the 
same purpose. 

In addition, we propose a simple way to calibrate the relationship between the ampli- 
tudes of the radial velocity and the beaming flux variation. Since the satellite motion is 
known, and is linked with the motion of the Earth, we already have a well known radial- 
velocity signal in the data for all stars. C0R0T, for example, will observe dense fields around 
the ecliptic, continuously for almost half a year. Thus, the amplitude of this heliocentric 
velocity signal will be close to 60kms _1 . 

The beaming photometric signal associated with the motion of the telescope will affect 
all stars, binary and single alike. Measuring this signal can serve to calibrate a with the 
radial-velocity amplitude, which in turn can be used to interpret the beaming signal of the 
stellar orbital velocity, if it is a binary. For single stars a is actually a piece of spectral infor- 
mation that reveals the location of the passband along the blackbody radiation curve, and 
thus provides an estimate of the stellar effective temperature. Figure [3] shows the expected 
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photometric variability amplitude in V for different temperatures, due to heliocentric mo- 
tion alone, assuming blackbody radiation law, and a heliocentric radial-velocity amplitude 
of 60kms _1 . 

Figure H] shows a simulated lightcurve that demonstrates the type of signal we expect 
to detect for a binary. The lightcurve includes the beaming variability of a 10-day period 
GO-KO binary star, together with the beaming variability related to the heliocentric motion. 
The associated radial velocity amplitudes are Ki = 52kms~ 1 and K 2 = GOkms™ 1 . The 
noise included is only a 10 -4 white noise. Note that in the case of the heliocentric motion 
the contributions of the two components are in phase, and therefore do not cancel out. 

Since the signal related to the observatory motion is common to all the observed objects, 
it will appear as a systematic effect, and may be mistakenly removed as such. Algorithms 



such as SysRem (ITamuz. Mazeh &: Zuckerl 120051 ) should identify such effects and also the 



individual response of each object to the same effect, and are thus ideally suited to provide 
the required information for its correct analysis. 

In the lightcurve we show in Figure H] we have neglected stellar microvariability. When 
actual lightcurves are analyzed, this variability should be properly accounted for. For 
older solar-type stars of low chromospheric activity, we may use the solar microvariabil- 
ity as the only available example . Although the solar rotational period is about 26 days, 



Aigrain. Favata &: Gilmord (12004 ) show that due to the short lifetimes of the spots and fac- 
ulae, there is no clear periodic signal in this period. Instead, most of the signal is in higher 
frequencies, corresponding to periods shorter than 10 days. For non-solar type stars, vari- 
ability due to chromospheric activity may be larger than solar, and care should be taken in 
separating the beaming effects from the variability effects. 

Efforts are currently underway to characterize the microvariability that we expect to 



observe with the high-precision p hotometric satellites (e.g.. lAigrain. Favata Gilmordl2004 



Lanza et al.l 120061 ; iLudwigi 120061 ). in order to facilitate the detection of planetary transits. 
However, the stellar microvariability is not expected to be strictly periodic, and therefore 
it should be possible to single out the beaming effect, specially because the shape of the 
beaming modulation is known and depends only on a few parameters. Spectroscopic infor- 
mation obtained as part of the follow-up observations should also help to further study the 
chromospheric activity of the binary candidates. 

In order to give an order-of-magnitude estimate of the expected number of beaming bi- 
naries that will be detected by, e.g., CoRoT, we estimate that 10% of the observed late-types 
stars have been discovered to be spectroscopic binaries with a threshold of K\ > 3kms _1 
and periods of less than a year or so. This estimate is based on results of the seminal works 
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of iDuqennoy &: Mayorl (Il99ll ) and lLatham et al.l (120021 ) . During the lifetime of the mission, 
CoRoT is expected to monitor about 60000 stars. Thus, we roughly expect 6000 of them to 
be detectable as beaming binaries with 10/zmag variability or higher. Since this estimate ap- 
plies to binaries with late-type primaries, and accounting for the fact that beaming is biased 
against equal-mass binaries, we can somewhat scale this number downwards, and estimate 
CoRoT to yield at least 1000 beaming binary stars. The fact that the CoRoT sample of 
beaming binaries will be discovered by a systematic, magnitude-limited survey, will augment 
significantly the statistical knowledge of binaries. 



4. Conclusion 



As we have shown in Section [21 we expect the multitude of very precise light curves 
of CoRoT and Kepler to yield hundreds of new binaries through their periodic beaming 
variability. Thus a new observational category will emerge - - beaming binaries. In all 
types of binaries, the discovery and the analysis of the binary motion strongly depend on 
the timing and the number of the measurements. For the satellite photometric data we 
expect continuous radial-velocity data that will yield the spectroscopic orbital elements, the 
period and eccentricity in particular. Once the spectral index a is known, the radial-velocity 
amplitude can be derived as well. 

Without the beaming binaries, CoRoT and Kepler are supposed to find most of the 
eclipsing binaries. However, binaries with periods longer than a few days need very fortuitous 
geometrical si tuations to present eclipses, an d are therefore rare in the data of photometric 



Ttan,., ft Nortj m. The beaming binaries can be detected up to 



surveys (e.g. 

a hundred days and more. Therefore, the new class of binaries will extend our detailed 
knowledge of the statistical characteristics of binaries by an order of magnitude, specially 
because these binaries will emerge in the context of a well defined, complete, magnitude- 
limited, photometric search. This can shed light on the distribution of orbital period (e.g., 
Duqennoy fc Mayorl Il99ll ; iMazeh. Tamuz fc North! 120061) and the relation between orbital 



period and eccentricity (e.g.. lHalbwachs et al. 



20031 ). 
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Table 1. Stellar parameters for simulation 



parameter 


F0 


GO 


K0 


Reference 


Mass [M Q ] 


1.6 


1.05 


0.79 


1 


Radius [Rq] 


1.5 


1.1 


0.85 


1 


Effective temperature [K] 


7300 


5940 


5150 


1 


Gravity darkening coefficient in V 


0.9 


0.4 


0.4 


2 



References. — (1) I Cox J2000h (2) iMorrisj (|l985h 



- 11 - 



Table 2. The three periodic photometric effects for sample binary configurations 



P = 10 days P = 100 days 

Primary Secondary ellipsoidal reflection beaming ellipsoidal reflection beaming 



F0 GO 3.9 -10~ 4 4.8 • 10~ 4 6.4 • 10" 4 3.9 • lO" 6 2.2 • 10~ 5 2.9 • 10" 4 

F0 K0 3.4 -10" 4 4.1 -KT 4 8.3 • KT 4 3.4 • KT 6 1.9 • 10~ 5 3.8 • KT 4 

GO K0 1.9 -10- 4 2.1-10- 4 6.6 -10- 4 1.9 • 10 6 9.6 • lO" 6 3.1 • KT 4 



- 12 - 



10 




Period [days] 



Fig. 1. — The three periodic photometric effects for an F0-K0 binary star in a range of 
periods. 
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Mass Ratio 

Fig. 2. — The three periodic photometric effects for a GO primary in a range of mass ratios. 




Fig. 3. — The amplitude of the expected beaming variability of a black-body radiation 
induced by the heliocentric motion of CoRoT, in the V band, assuming the star location is 
close to the ecliptic and the radial velocity amplitude is 60kms _1 . 
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Fig. 4. — A simulated CoRoT lightcurve of a 10-day period G0-K0 binary star. The beaming 
effects related to the binary orbit and the heliocentric motion are easily noticed. The dashed 
line represents the pure heliocentric motion beaming signal. 



